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ABSTRACT 

The luminosity function of quasars and its evolution are discussed, 
based on comparison of available data on optically selected quasars and 
quasars found in radio catalogs . It is assumed that the redshift of quasars 
is cosmological and the results are expressed in the framework of the A = 0, 
q = 1 cosmological model . The predictions of various density evolution laws 

are compared with observations of an optically selected sample of quasars 

J 

and quasar samples from radio catalogs. The differences between the optical 
luminosity functions, the redshift distributions and the radio to optical 
luminosity ratios of optically selected quasars and radio quasars rule out 
luminosity functions where there is complete absence of correlation between 
radio and optical luminosities. These differences also imply that Schmidt's 
( 1970 ) luminosity function, where there exists a statistical correlation 
between radio and optical luminosities, although may be correct for high 
redshift objects, disagrees with observation at low redshifts. These 
differences can be accounted for by postulating existence of two classes 
(i and II) of objects. The class I non-evolving objects dominate the opti- 
cally selected samples at low redshifts, have low optical luminosities of 
22 0 

about 10 w/Hz (at 2500 A) and are absent from radio samples because of their 
relative low radio to optical luminosity ratios. The class II objects are 
the quasars found in radio catalogs, show strong evolution, have optical 

23 24 

luminosities of 10 to 10 w/Hz and dominate the optically selected samples 

at large redshifts. Possible relations between the two classes are also 

discussed, and tests for the validity of this picture are presented. 

* 

Alfred P. Sloan Foundation Fellow 



I . Introduction 


If the redshift of the quasar is cosmological then their distribution 
in space-time indicates that the universe is not in a steady state and that 
the quasars show strong evolution with cosmic time in conventional Friedmann 
models. This also means that the quasars cannot be used for the determina- 
tion of the parameters (such as deceleration parameter, etc.) of the model 
for our universe. The subject of this paper is the evolution of the lumi- 
nosity function of the quasars in the framework of a given cosmological 
model with the assumption that their redshift is cosmological. Unless 
otherwise stated, we shall assume the cosmological model with zero cosmo- 
logical constant and with a deceleration parameter q = 1. In Section II, 
we derive equations relating the luminosity function and its evolution to 
the observations of quasistellar radio sources (quasars observed first at 
radio frequencies and identified with optical objects subsequently) and 
optically selected quasars (or the so-called radio quiet quasars). The 
luminosity function and its evolution can be described by a function 
HL,F,z) which gives the frequency distribution of the optical luminosity 
L and radio luminosity F of objects at different cosmological epochs 
or redshifts z. The results of comparison with observation are discussed 
in Section III where we compare various density evolution laws, and in 
Sections IV and V where we discuss various possible luminosity functions 
and their extrapolation to low luminosities. In Section VI we present a 
modified V/V^ analysis for optically selected quasars. The results are 
summarized in Section VII. 
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1 1 . Evolution of the Luminosity Function 


The luminosity function and its evolution can be represented by a 
single function Y(L,F,z), which gives the frequency distribution of ob- 
jects in the three dimensional space L,F, and z. [Following Schmidt (1968) 

L will refer to the luminosity at 2500 A (Schmidt's F(2500)> and F will 
refer to the luminosity at 500 MHz (Schmidt’s F(500)).] This function gives 
an average statistical description of the variation of the optical and radio 
luminosities of the objects throughout their lifetime, their birth and ex- 
tinction rates and the changes of these with cosmological epoch. It is 
hoped that the knowledge of the luminosity function ¥ will give us some 
understanding of the physical processes which govern the behavior of these 
objects. The function Y itself is deduced from observations of the red- 
shift, optical intensity I, and radio flux density S of the objects. 

In order to account for some selection effects, one needs this data for a 
sample of sources complete to a given optical intensity I q and radio flux 
density S q . Given this data one can calculate the optical and radio lumi- 
nosities of each source in the framework of any cosmological model, 


F = Sg r (z) « Sz 2 , L = Ig Q (z) Iz 2 . 


( 1 ) 


Note that because we are interested in luminosities at a given frequency 

the functions g and g will be different due to the difference in 
r o 


"*C^r 


"Ob 


the optical and radio spectra of the sources. If S a v '~ i ', lav , 
then for our cosmological model 


g x (z) = 4 jt 




z 2 (l+z) 1 ax 


( 2 ) 


Because (a r ) ~ (a ) ~ 0.7, the K-correction term (1+z) 1 is small 

and the same at radio and optical frequencies. We shall neglect this 

2 

term in our discussion and set 4jt(c/H Q ) =1 as indicated by equation 
(1) . However, the numerical results are obtained including the K-correc- 
tion terms and with (c/H ) = 3000 Mpc . Note also that I is same 

o o 

as F (2500 ) of Schmidt (1968). 
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Then from Nj o> s 0 (F*L), the observed number of sources in the sample with 
optical and radio luminosities L to L + dL and F to F + dF we have 


f m 

J Y(L,F,z) V' dz = Nj g (F,L) , v ra = v ( z m ) ( 3 ) 

O 0*0 


where V(z) is the co-moving volume up to a redshift z [for our model V ; 
2 

dV z 

-r- oc -= 1 and z is the maximum redshift within which the objects 

dz . N 3 n; — — m 

(l+z) J sll+2z 

of luminosities L and F could be observed given the limits I ,S . 

o o 


z 

m 



z 

m 



j 



radio limited 


o 



optically limited . 


o 


(*) 


We can also write 

dN s (L,F) 

* y ^ O r O 

Y(L,F,z ) « 2__2 x (1+z ) 5 sll+2z /z ( 5 ) 

' ’ J m 7 dz v nr v m m v ' 

m 

so that by comparing the number of sources in samples limited to differ- 
ent I ,S we can calculate Y. 
o 7 o 

Unfortunately, with the limited data available at the present time one 
cannot deduce a complete description of the distribution in L,F and z. 

For example, the 3CR sample of Schmidt ( 1967 ) and the 4C sample of Lynds 
and Wills (I 971 ), each contain about 30 objects. Subdividing these 
samples into various luminosity classes and various limiting observed 
intensities (or maximum redshifts, z m ' s ) give only few objects in each 

subgroup and reduce the statistical significance of the results. Further- 
more, since low luminosity objects can be seen only when they are nearby we 
cannot learn about the behavior of these objects at large redshifts. The 
high luminosity objects, on the other hand, are less numerous and we must 
sample a larger volume of space to observe them. The probability of finding 
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these objects at small redshifts is less than the probability of seeing 
them at large redshift so that one cannot learn about the behavior of ¥ 


for the high luminosity objects at low redshifts. Because of these limita- 
tions it is necessary to compare some average integral property of the 
function Y with observations. The one moment of Y used in recent years 
for comparison with observation, which renders to simple interpretation in 
special cases, is the average value of V/V m [suggested first by Kafka 
(1967) and applied to 3 CR sources by Schmidt (1968) and by Lynds and Wills 
( 1972 ) to 4 C sources] . In the Appendix we show the expression for 
for quasistellar radio sources which have two limiting fluxes, radio and 
optical. However, the essential feature of this test can be seen more clearly 
if we consider a sample of sources with only one limiting flux. This will 
be the case for a sample of optically selected quasars where only the opti- 
cal limit must be taken into account . If we integrate over all radio lumino- 
sities to obtain the optical luminosity function 


*op t (L ' z) = J Y(L,F,z)dF , 

o 


(6) 


then the average value of V/V m is given by 


< x ) 3 ( V/V m) ‘ ~ 


00 ^ 

/ / 


m 


W opt ( L ,a)v'dz 


eo V 

J dL J 

o o 


( 7 ) 


Y opt( l ’ 2i)v ' d ' 


It is clear that if the luminosity function does not evolve (y (L,z) 
independent of z) then = 0.5. In general, if the observed ^x^ = 0.5, 

evolution cannot be ruled out. But if ^x^ ^ 0.5 then there must be some 
evolution of the luminosity function . In order to proceed further and 
determine the possible evolutionary laws we must make some simplifying 
assumption about the luminosity function. The various ways that the lumino- 
sity function can be written and the interpretation of these ways have 
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been discussed elsewhere (Lynds and Petrosian 1972) . We shall limit our 
discussion here to the simplest type of evolution, namely the case of pure 
density evolution where it is assumed that the shape (dependence on L and F) 
of the luminosity function is epoch independent, so that one can write for 
the general case 


Y(L,F,z) = Y(L,F)p'( z) 



L,F)dLdF 


1 , 


( 8 ) 


and for the case of one limiting flux (equation (6)) 

' f opt( L < a ) - ,f op t < L >e( z ) ’ f f (L)dL , 1 . 

J opt x ' 


(9) 


With this assumption, expressions (3) and ( 5 ) for the luminosity function 
are simplified as 

r V m(L,F) 


/f 1 

Y(L,F) = N s (F,L) / j p(z)v'dz 

o’ o / IJ J 


( 10 ) 


and for the case of one limiting flux 


t.ptW ” N 


W[/ 


m 


p(z) V 'dz 


( 11 ) 


where Nj q (l) is the number of sources with luminosity L to L+dL in an optically 
selected sample of sources with I > I Q . Substitution of equations (9) and (11) 
into equation ( 7 ) gives 


N 


00 o 

<*> = J N^(L)x(L)dL = ^2 X ( L i) > 


( 12 ) 


o o 


i=l 


where 


J m 

Vp(z) v'dz 

x(L) = 


v m J ra p( z ) v ' dz 

o 


(13) 
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00 



/ 


Nj (L)dL 
o 


is the total number of sources in a sample complete to a 


o 

flux density I , and z = vjL/I . 

J o m v o 

From observed values of x(L) and ^x^ and equations (12) and (13) 
one can determine the law of the evolution of density p(z). Then with 
the help of this law and equation (10) one can calculate the luminosity 
function. 
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III. Evolution 


On Table 1 we summarize the observed values of ( V/V ) for the 3CR 

sample (Schmidt 1968), the 4 C sample of Lynds and Wills (1972) and a sample 

of optically selected objects (Lynds et al (1972)). (The V/V values for 

m 

the latter sample were obtained by a slightly different method. This method 

is described in Section VI.) The first number gives the value of (x) for 

the whole sample. The second and third numbers give the values of ( x) for 

the high and low (or luminosity) objects, respectively. The number of 

sources, N, and the expected standard deviation a = — r == • are also given 

V12N 

for each sample. No useful information can be obtained from further subdi- 
vision of the samples because of large values of the o's. 

It is clear that with such meager data one cannot give a unique descrip- 
tion of the function p(z). As a result a few single parameter functions 
have been considered. We shall consider here the following functional forms 
for p(z). 



p(z) 

= p o (l + z) ra 

(14) 


p(z) 

a (V(z)) a 

(15) 

p(z) 

= P o 

exp (px(z)) 

(16a) 

p(z) 

= P o 

exp {p’z/(l+z)j 

(16b) 


where t(z) is the look back time in units of Hubble age h 1 . For the 
cosmological model under consideration here 

t(z) = 2 tan 1 n/i+2z + Vl+2z/(l+z) - (l+rc/2) . ( 17 ) 

The redshift dependences of the various evolution laws are shown in Figure 

1. Figure 2 shows the variation of x(z ) with z (which is equivalent 

m m 

to the variation of x(L) with L) for the above density evolution laws 
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(equation 13). An evolution law nearly identical to a power law was first 
suggested for quasars by Schmidt. From theoretical considerations a power 
law in (1+z) appears reasonable on the grounds that most physical quanti- 
ties (such as average mass density etc.) in the universe, which could pos- 
sibly effect the density of objects vary with some power of (1+z) (see, 
for example, Rees and Setti, 1968). However, this law clearly cannot be 
extended to large redshifts because it predicts more objects with redshift 
>2.5 than has been observed (Petrosian, 1970) and because it would result 
in an isotropic background intensity at meter wavelengths much larger than 
observed. Furthermore, as it was pointed out by Lynds and Wills (1971), 
for the 3CR and 4C samples the observed values of ( x) for high and low 
luminosity objects are equal, while the expected values of ( x) from a 

power law density evolution for samples with the same z (or luminosity) 

m 

distribution as the 3CR and 4C samples are different (cf. Table 1). 

Although for each sample alone the differences between the expected 
and observed values of ( x) are not statistically significant, the fact 
that both samples show similar behavior was considered significant by Lynds 
and Wills (1972) who suggested the density evolution law of equation (15) 
for which (as is evident from equation (13)) x(L) is independent of L and 


<*> 


a + 1 

a + 2 ' 


(18) 


This is also true for the general luminosity function (see Appendix 1). 

The unsatisfactory features of this evolution law are: (i) it implies zero 

density of objects at the present epoch, and (ii) it is difficult to find a 
physical relation between the co-moving volume and the density of sources. 

The objections to the above evolution laws could be overcome by slight 
modifications. For example, modifying equation (14) so that the exponent 
m is a decreasing function of redshift will overcome the difficulties 
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encountered with that evolution law. Similarly if equation (15) is modified 

OL 

to read p(z) = a + b V , the present density of objects in the universe 
becomes finite. 

Exponential evolution laws similar to equation (l 6) have been suggested 
on various occasions. This law was used first by Rowan-Robinson ( 1967 )* and 
later by Doroshkevich et al . (1970) in connection with the evolution of radio 
sources. In an earlier paper ( 1970 ) I have shown that an exponential 
law of equation (l 6 b) fits the redshift distribution of quasars with- 
out requiring an additional cutoff at larger redshifts. As is evident from 
Figures (l) and (2), the exponential laws predict slow variation of 
x(l) with L. Table 1 shows the expected values of (x) for the evolution 
law of equation (l 6 a) which are in better agreement with the observed values 
than the power law. Lynden-Bell (1972) has arrived at essentially the same 
conclusion using a different analysis of the 3 CR data. Thus, the exponen- 
tial laws overcome the objections raised against the previous two laws. 
Furthermore, an exponential evolution law could come about if the processes 
of the physical evolution of the sources (i.e. variation of their luminosity 
with their proper time) are statistical in nature as in nuclear radioacti- 
vity. Such a model has been discussed by Bahcall ( 1972 ) . 

Let us now consider the optically selected sample of Lynds et al . (I 972 ). 
As is evident from Table 1, for this sample, the values of (x^ for the 
high and low z^ objects are different'. This is in good agreement with 
the expected values from a power law density evolution but shows statisti- 
cally significant deviation from the expected values of the exponential laws 
or the law of equation ( 15 ) . This evidence clearly reduces the significance 
of the disagreement between the observed behavior of the 3CR and ItC samples 
and the expected behavior of the power law density evolution. However, if 
the 3CR and 4C data are taken as evidence for the exponential law of 
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evolution (or that of equation (15)), then we must modify these laws at low 
redshifts. Since the value of (x) for the low z m objects in the optically 
selected sample is not significantly different from 0.5 we may modify equa- 
tion (15) or (l6) by adding a constant term on the right hand side. For 
example, if we modify equation (15) to read 


P(V) 



then from equation ( 13 ) we have 


x 



1 + (V /V ) a (l+o/2) 
Qfl-1 o nr v ' 

af2 1 + (v /v ) a (id-a:) 

x o nr x ' 



V 

m 


« V(z ) 
v o 


V » V ( z ) 
m o' 


( 19 ) 


which by proper choice of z^ (z q « 0.6) could give the observed differ- 
ence between the values of (x^ for high and low z^ objects. Note 
that this also overcomes one of the objections raised against this evolu- 
tion law earlier. Similarly equation (l6a) can be modified as 


P( z ) = Pj + PI / T(Z) • (20) 

Note that these modifications imply the existence of two classes of 
quasars. One class (class I) which shows no evolution and consequently is 
observed primarily at low redshifts. A second class (class II) which at 
the present epoch are less numerous than the first class but which have 
evolved such that at earlier epochs they were more numerous. Note also that 
the ten out of the eleven objects with low values of z m in the optically 
selected sample have measured redshifts z < O.k while the average red- 
shift of the other 12 objects is about 1.4. As we shall see in Section 
V there appears to be other indications of differences between high redshift 
( z ~ 1 to 2) and low redshift (z < 0.6) objects. 
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In summary, if we assume that quasars have undergone density evolution 

then the combined V/V analysis of 3CR, 4c and optically selected quasar 

m 

data shows that a power law density evolution, although not the best fitting 
law, nevertheless is acceptable within the statistical uncertainty of the 
data. On the other hand, if we accept the 3CR and data as evidence for 
exponential type evolution, then the data on radio quiet objects indicates 
the existence of two classes of quasars. 
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IV. The Luminosity Function 


Assuming density evolution then equation (10) and (ll) and the density 
evolution law obtained from the V/V m analysis give the luminosity func- 
tion. If the assumption of density evolution is correct the above lumino- 
sity function would be valid for all epochs, otherwise it only would be 
an average luminosity function for the range of the observed redshifts. 

For samples consisting of small number of sources these equations must 
be integrated over a considerable range of luminosities so that it in- 
cludes contribution from a few objects. Thus in practice one calculates 


V^F) = 


i L X A ~ F X A m ■ ^ 

dL'J dF'Y(L'lO = y] 

*F i 


/ 


( 21 ) 


mi 


p (z) dV 


o 

where the summation includes contribution from all objects with observed 
luminosities in the indicated ranges. This procedure has been applied to the 
3CR and the 4c sample by Lynds and Wills (1972) who present their results 
(for A = 2.5) in tabular form and plot $ 1 (F,L) versus F for various 
values of the optical luminosity and versus L for various values of the 
radio luminosity. It appears from these plots that the radio luminosity 
function can be fitted to a power law with a power exponent which is 

approximately the same for all optical luminosities. The optical lumino- 

23 8 

sity functions are flat (with a steep cutoff for L ^ 10‘ w/Hz) and 
show larger dispersion in the power exponent. If the above exponents were 
indeed independent of radio and optical luminosities we could treat L 
and F as stochastically independent parameters and write 


T ( L ’ F > - ,f rad( F > • < 22 > 

In order to reduce the scatter in the data we have plotted in Figure 3 , the 
cumulative luminosity function 
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(23) 


*(L,F) 



Y(L , ,F 'JdL ' 


versus F and L for various values of L and F respectively (because of 
steepness of the luminosity function the cumulative luminosity function 
is nearly identical to the function $ in equation 21). As can be seen 
from these graphs identical curves could be fitted to all these points 
(although the curves plotted are not the best fitting curves) implying 
that equation (16) is approximately valid. The curves shown in Figure ( 3 ) 
can be expressed analytically as 




ft) 


- 1.3 


-e 


F > F =10 
o 


25.9 


w/Hz 


[l+(L/L o ) ] 


( 24 ) 


L>10 1,6 L , L =10 23,8 w/Hz 

o o 


which approximately corresponds to a differential luminosity function of the 


form 


Y(L,F) oc F 2> 3 L 1 '® , 

[l+(L/L o ) n ] 


/ F > F 

II > 10 1,6 L 


( 25 ) 


The curves on Figure (3b) are for e = 0.3. The data allow values of e 

ranging from 0 to 0.5. The term [l+(L/L Q ) n ] with n > 1 describes the 

24 

apparent cutoff for L > 10 w/Hz. No such term seems to be necessary 
for high radio luminosities. 


Of course there are no compelling reasons to write the luminosity 
function as in equation (22). For a small sample of sources, it is possible 
to find other uncorrelated pairs of variables X = f 1 (L,F), Y = f (l,F), 
in which case one can repeat the above analysis (equation ( 21 ) onward) by 
changing L and F to X and Y. For example, Schmidt (I 97 O) has suggested 
that the radio and optical luminosities are statistically correlated so that 


one must write the luminosity function in terms of the ratio R = F/L and L 
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in which case instead of equations (22) and ( 23 ) we have a new luminosity 
function 


Y'^R) = ^ pt (L)^(R) ' 


(26) 


. CO 00 


f'(L,R) = 


f dL'f dR / Y / (L / ,R / )= 

•fe i j* mi p 


( 27 ) 


(V)dV 


As before, equation ( 27 ) expresses the new luminosity function in terms 
of the observed quantities, so that it can also be derived from Table 3 of 
Lynds and Wills (1972). The results of this derivation are shown on Figure 
4 where 4>' is plotted for various values of R and L versus L and R respectively. 
The two representations (equations 22 and 26 ) are related as 


¥'(L,R)dLdR = Y(L,F)dLdF , 

so that the new cumulative luminosity function is given by 


(28) 


* , ( L,S ) ~ E_1 ' 3 / dL ' L ,a.3nr i 1 


1 + 

(¥-\ n 


i L o) J 


( 29 ) 


-1 .3 

For a given L, <&’ <x R . Curves representing this relation are shown in 

Figure (4a). These curves agree with the observed points to the same accur- 

, -1.3-e 

acy as the curves in Figure (3a). Similarly, for a given R, 0 (R, L) <x L 

for L « L (dashed lines in Figure (4b) plotted for e = 0.2) and $'(L,R) <x 
o 

“"H”! 3_ e 

L " for L » L . These relations and equation (29) can be approximated as 

o 


$ / (r,l) oc r- 1 -^ 1-3-e 


1+(L/L q ) 


L S 10 _1 ° 6 L 

3.0 


0 


( 30 ) 


R > 10- 

Lynden-Bell (1972) using a different method for the 30R sample alone arrives 


at a similar luminosity function. 



Equations (24) and (30) are two different representations of the data 


on the quasistellar radio sources. They agree with observations to the 

same accuracy over the specified ranges of the luminosities. It is clear 

that these functions cannot be extended to very low luminosities. There 

must be a flattening of the optical luminosity functions (in both cases) 

for low values of L, and there must be a break in the power law for low 

values of F (in equation (24)) and for low values of R (in equation 

(30)). The position of this break is the distinguishing feature of these 

luminosity functions. If for all values of L the break occurs at the 

same value of F (or R) the correct representation of the luminosity 

function is given by equation (24) (or by equation (30)). As is evident 

from Figures (3) and (4) one cannot distinguish between the two luminosity 

functions based on quasistellar radio source data (3CR and 4C samples) 

alone. For example, these luminosity functions predict nearly identical 

S 

redshift distributions for all samples with a ratio of comparable to 

4 ° 4 

that of the 3CR (S /I = 9 X 10 ) and the 4C (S /I = 5.6 X 10 ) sam- 

o o o o 

4 

pies. Only for samples with the ratio S /I « 10 , in particular for 

o o 

any sample of optically selected quasars (S /I = 0), do the predictions 

o o 

of the two luminosity functions differ substantially. We compare the ob- 
servations of optically selected quasars with the prediction of the two 
luminosity functions in the next section. 
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V. Extent ion of the Luminosity Function to Low Luminosities 


and Optically Selected Quasars 


In this section we compare the observed redshift distribution of quasars 
with the prediction of the two luminosity functions discussed in the pre- 
vious section. For the case of stochastic independence of radio and optical 
luminosities (equation (22)) we adopt the cumulative luminosity function 


i . 

opt 


i(L,F) = » (LKf/F }- 1 ' 3 , 


opt 


< L ) =/( L/L mi»>" E 


1+10 


-2n 


_ 1+ ( L/L q ) 


F > F = 10 2 5-9 w/Hz 

O 


L > L . = 10 \ = 10 21 ’ 8 w/Hz 

mm o 

L ’ n (31) 


where we have not specified the form of the radio luminosity function for 

25 . 9 

sources with luminosities (at 500 MHz) less than 10 * w/Hz and have as- 

sumed absence of sources with optical luminosities (at 2500 A) less than 
21 . 8 

10 ‘ w/Hz. These numbers correspond to the lowest luminosities among 

all objects with known redshifts. The frequency distribution of redshifts 
of objects brighter than optical flux density I q and radio flux density 
S q is given by 

N(z) = V ' P(z) $ opt ( V 2 ) * rad (V 2 ) • (32) 

In Figure (5a)we present schematically the variation of this distribution 
with redshift. The top curve presents ' p(z) which is the redshift dis- 
tribution of all observable sources (I = S = 0). For optically selected 

objects (s = 0. $ , = l) we obtain the lower solid curve, which begins 

1/2 

to deviate from the upper curve at the redshift z = (L . /I ) .At low 

opt v min o' 
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values of z the deviation is small and is a slowly varying function of 

z because of the flatness of the optical luminosity function. The rapid 

1/2 

cutoff at large z [(z > (L /I Q ) ], is due to the steepening of the 
luminosity function for L S L q . 

For a sample of quasistellar radio sources (S > 0) we must also 

' o 

include the radio selection effects which is given by 

$ (S z 2 ) = (z/z ) 2 *^ for z > z = (F /S ) //2 (33) 

if there are no sources with radio luminosity F < F . The dashed lines on 

Figure (5a) give the expected redshift distribution for two such samples, the 

top line for the case S /I < F /L . and the lower line for S /I > F /L 

o o o min o o o min 

4 \ 4 

as 10 . As mentioned in the previous section, in practical cases S /I q > 10 

so that the lower case is more commonly applicable. 

As it has been pointed out earlier (Petrosian, 197°) these redshift 
distribution curves agree with the observed distribution of quasistellar 
radio sources (3CR, and all known objects) but not with the redshift distri- 
bution of all known optically selected quasars, which has a double peaked 
distribution. As we shall see below this disagreement cannot be entirely 

accounted for by observational selection effects, so that this luminosity 

4 

function cannot be applied to samples with S /I < 10 , unless it is modi- 

o o 7 

fied, for example, as in equation (40) below. 

Furthermore, analysis of optically selected samples of quasars (Braccesi, 
1972, Schmidt, 1970) show that the optical luminosity function is much steeper 
than that given by equation (31) (compare Figures 3b and 6), which also speaks 
against the validity of this luminosity function for low radio luminosities. 
The observed steep luminosity function for the optically selected quasar would 
come about if the optical luminosity function becomes steeper for low values 
of radio luminosities. This means that e in equation (31) is a function of 
the radio luminosity and it increases with decreasing F. There is, however 

17 



no evidence for any systematic steepening of the optical luminosity function 

■with decreasing F for F > F . To determine whether rapid steepening of the 

optical luminosity function takes place for F < F we must wait for data 

o 

on samples with more objects and extending to lower values of S 

o * 

A steep optical luminosity function would be observed if, as suggested 
by Schmidt (I 97 O), the radio and optical luminosities were correlated so 
that the form of the optical luminosity function is independent of the 
value of R = F/L. Then we have equations (28) and (30) or 


§'(L,R) = S 0 ' pt (L) G(R) 


, L < L . = 10 _2 L 

min o 


$' .(L) = 

opt v ' 


(L/L ) 
' min' 


-1.3-e 


1+10 


- 2 n 


1 +(L/L ) J 


, L > L 


min 


G(R) = K(R/R ) 1,3 , R > R = 10 3 ' 2 


G(R . ) = 1 
' min' 


(34) 


where we have not specified the value of R . or the form of the function 

mm 

G(R) for R > R > R . . 

o mm 

In this case the redshift distribution (for z > z . ) of objects 

v opt ’ 

brighter than I and S is: 

00 

i) 0 < S < I R . 

o o min 


N(z) = p(z) v' 


, z < z = (L . /I ) 
’ opt v mm o' 


1/2 


$ / (I z ) , z > z 
opt v o ' ’ opt 


ii) S > I R . 
o o mm 


, z<z =(L.R./S) 
’ r min mm o' 


(35) 

1/2 


N(z) = p(z)v' 


/S z 


^optlR . )[l+X(z) , z > 

\ min/L J 
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where 


X(z) = 


(k) ^ 


S Z 2 /R . 

° " ln /S z 2 \ 

’foV L > G (-H 

I z 2 
o 


dL 


(35)' 


is a function of order unity changing slowly from zero at z = to a 
constant value for z S z As s ^ own schematically in Figure (5b) the 

redshift distribution for the quasistellar radio sources (i.e. S /I > R . ) 
expected from equation (3*0 is identical to that expected from equation 
(32). However, for S q /I o ^ (i- e « for optically selected quasars) 

the redshift distribution of equation (35) i s nearly identical to the 
redshift distribution of quasistellar radio sources, because in this case - 
unlike in equation ( 31 ) - the optical luminosity function ^Q pt (b) is as 
steep as the radio luminosity function ^ rac j(^)* * n the luminosity 

function of equation (3*0 was proposed by Schmidt ( 1970 ) when he observed 
a rough agreement between the redshift distribution of a sample of 18 magni- 
tude optically selected quasars and the 3CR quasars.* However, as we shall 
discuss below, data from other optically selected samples of quasars does 

not agree with some of the predictions of this luminosity function. In ad- 

dition, this luminosity function encounters the following difficulty. 

Comparison of observed and expected values of R : If equation (34) is 

valid for all radio luminosities (or for all values of R), then for any 

sample of optically selected quasars the average value of the ratio of radio 

to optical intensities (or luminosities) is given by 

(36) 


(!) - (!) - ( r ) - f RG '< E > dB > 

R . 
min 

where G^R) = -dG(R)/dR. There has been various attempts of radio 


* 

The independence of redshift distribution from the radio flux density limit 
may also be considered as evidence for non-cosmological nature of the red- 
shifts . 
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observations of optically selected quasars. The most recent attempt by 
Wardle and Miley ( 1971 ) indicates positive detection of radio emission (at 
2695 MHz) only from one object. For all other objects of known redshift 
they failed to detect any radio emission greater than a few times 10 ^w/(m^Hz) . 
Using the upper limit for each source from these data and assuming an average 
radio spectral index C£ r ~ 0.7 we find that 


< R >, 


observed 


< 200 


(37) 


The expected average value of R from equation (3^-) will depend on 

q p 

the behavior of the function G(R) for R <, R q = 10 0 * . Schmidt (1970), 
from comparison of 3CR and radio quiet quasars, gives G(R=R q ) = K ^ 0.05. 

It can be shown that the smallest expected value for ( R) is obtained if 
G(R) °c log R for R < R q , so that the expected value of ( R) for the 
luminosity function of equation (34) is 

< E >eq . 34 > 5°° • (3 8 ) 

Note, for example the extention of G(R) suggested by Schmidt (I 970 ) gives 

O 

<R> fa 10 . The sample of all known optically selected quasars includes 
many PHL objects, few of which are identified with 4C sources. Since the 
redshifts of all PHL objects are not known it is not clear how much the 
identified PHL objects contribute to the value of (r^ in equation (38). 

On the other hand, all the objects in the sample of Lynds et al . (I 972 ) 
were included in Wardle and Miley’ s observation program so that for the 22 
objects in this sample 


( R > observed < 120 

compared to the expected value of equation (38). The discrepancy between 
equation (39) (or equation (37)) and equation (38) may be due to the 
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difficulty in the radio observations of these objects. Wardle and Miley 
indicate that few of the objects are confused with strong radio sources. 

If it just happens that these objects are themselves strong emitters, the 
above discrepancy could be resolved. In a recent paper, Schmidt (1972) 
refers to observations by Katgert et.al. (unpublished) who have detected 
radio emission at 1410 MHz from 2, or possibly 3, objects at the level of 
0.01 flux units. If this data is confirmed, the above difficulty with the 
luminosity function of equation (34) will be alleviated. Otherwise this 
equation must be abandoned, partially or completely as a description of 
low radio luminosity quasars. 

There are various ways of extending the observed luminosity function 
of quasistellar radio sources to lower radio luminosities which could over- 
come such discrepancies. For example, if the observed redshift distribution 
of optically selected quasars and their steep luminosity function is due to 
the steepening of the optical luminosity function of equation (31) for F<F 

o 

rather than due to the existence of a correlation between optical and radio 
luminosities, then the luminosity function is given by equation (31) for 
F > F q and the optical luminosity function by 

$ ( L ' F min) = f dL ' f dF'Y(L'F*)ce l" 1 ’ 3 " 6 , (40) 

*L F . 

mm 

where F min is the smallest radio luminosity. In this case the expected 
values of (F/L) for optically selected quasars is given by 



(41) 


For F min ~ 10 * W//Hz this ratio w °uld be in agreement with the low ob- 

served value in equation (39). 
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On the other hand, if our earlier hypothesis that most of the low 


redshift objects belong to a different class is correct (see also the 
discussion below), then the above comparison should be made for objects 
with z > .6; i.e., if we assume that the luminosity function in equation 

(34) is valid only for the high redshift objects, then the value of K in 
equation (34) is reduced by a factor of about 2 which is the ratio of num- 
ber of all objects to number of objects with z < .6. Furthermore, the 
large redshift objects are optically weaker (on the average by a factor of 
2) than the low redshift objects. Inclusion of these two factors would 
eliminate the factor of 4 discrepancy between equations (38) and (39) . 

Comparison with observed redshift distribution : A more serious diffi- 

culty with the types of luminosity functions we have been discussing is their 
failure to explain the presence of relatively large number of low redshift 
objects among optically selected quasars. The observational data showing the 
relative abundance of low redshift objects (z < 0.6) in various samples of 

fV 

quasars are summarized in Table II. For the complete 3CR sample (Schmidt, 

1968) and the 4C sample (Lynds and Wills, 1972) the fraction, f, of objects 

with z < 0.6 is less than 30$> while for the optically selected sample of 

Lynds et al (1972), which is complete to an infrared magnitude i = 17.65, 

11 out of 22 objects have redshifts less than 0.6. The average value of 

visual minus infrared magnitudes for the Lynds et al (1972) sample is 

(v-i) = 0.65* which implies a visual limit of v as 18.3. From Figure 2 

of Braccesi et al (1970), we estimate the difference between ordinary visual 

magnitudes m and their magnitude v to be m -v as 0.05, so that the 
v v 

limiting visual magnitude m^ for this sample is similar to that of the 
3CR sample (m^ as 18.4). There is, however, one selection effect which 

* ' 

For definition of i and v magnitudes, cf . Braccesi et al, 1970. 
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could account for some of the difference between these samples. The average 
v-i colors of the objects with z < 0.6 in the Lynds et al sample is 0.2 
magnitude larger than for objects with z < 0.6. Whatever the cause of this, 

it will produce an over-representation of low redshift objects by a factor 

0 • 20 X 0 . 75 

of 10 ‘ 1.4, where 0.75 is the observed slope of the number magni- 

tude counts (Sandage and Luyten, 1969; Braccesi and Formiggini, 1969). Or, 
if we exclude objects with 17.45 < i < 17.65 from the list for z < 0.6, 
the fraction f for the optically selected sample is reduced to 9/20 as 
shown in Table II. Other differences in techniques of selection of objects 
for these two samples are small to account for the difference between the 
values of the fraction f . As can be seen from Figure 5 or from examination 
of equations (35) and (32), a larger value of f for optically selected 
quasars compared to radio quasars is contrary to the prediction of both lum- 
inosity functions. It clearly rules out the luminosity function of equation 
(31). The luminosity function of equation (34) predicts nearly identical 
redshift distribution for all samples of quasars limited to the same optical 
flux density I q independent of their radio flux density limits. In fact, 
contrary to observations, this luminosity function predicts (see Figure 5b) 
fewer low redshift objects (i.e., smaller values of f) for an optically 
selected samples compared to a radio sample. 

In order to test the significance of the differences between the observed 

2 

and expected behavior of the fraction f, we have carried out x analysis 
comparing the redshift distribution of the Lynds et al (1972) sample (after 
correction for the above mentioned selection effect) with the distribution 
of the 3CR and 4C samples. For the 4C sample we have included only objects 
with log I q > -30.00, which corresponds to m v < 18™ 4 (cf . Schmidt, 1968). 
The difference between 4C and 3CR samples is not significant. The difference 
between the fractions given on Table II for the optically selected sample and 
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the combined 3CR and 4C sample is larger but is not statistically significant. 
However, the dividing redshift z = 0.6 was chosen to divide the complete 
Lynds et al (1972) sample into two equal parts. This sample shows more pro- 
nounced concentration of sources around z ~ 0.3, so that dividing the sam- 
ple into two groups with z < 0.4 and z > 0.4 we find a probability of 
0.06 that the difference between the optically selected sample and the com- 
bined 3CR and 4C samples is due to chance . Furthermore, if the samples are 
divided into three groups, each group containing approximately equal numbers 
of optically selected objects, we find a probability of less than 0.01 that 
the difference in the redshift distribution is due to chance. As a test for 
the validity of equation (34) these results indicate even smaller probabili- 
ties since as mentioned above this equation predicts fewer low redshift ob- 
jects for an optically selected sample than for a radio sample. 

Equation (34) predicts identical redshift distributions for samples of 
sources with a given optical intensity. Therefore, a better test of this 
equation will be the comparison between the redshift distribution of samples 
of sources within a narrow range of magnitudes. Since most of the sources 
in the observed samples have magnitudes around 18 magnitude, we consider now 

o 

only 18 magnitude objects, i.e., sources with -30.00 < log I < -29.6 w/(Hzm ). 
In this case we can also include in our test Schmidt's (I 97 O) optically 
selected sample. 

The individual redshifts of the objects in Schmidt's (1970) optically 
selected sample are not published. For the 18 magnitude objects in this 
sample the fraction f = 5/19. There is one 18 magnitude object whose red- 
shift is not known and there are three 17 magnitude (log I > -29.6) objects. 
Thus, as shown in Table II for the complete sample f could be as large as 
9/23, a number very similar to that of Lynds et al (1972) sample. On Table 
II we also list the fraction f for the 18 magnitude objects in all 4 samples. 
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For the Lynds et al sample we give the fraction f obtained by including 

only objects with 16.45 < i < 17.45 for z < 0.6 and objects with 

16.65 < i < 17.65 for z >0.6. This accounts for the selection effect 

2 

discussed above. X analysis shows that the difference between the two 
optically selected samples is not statistically significant (the probabil- 
ity that this difference is due to chance is greater than 0.1). However, 
the probability that the difference between the combined optically selected 
sample and the combined samples from radio surveys (3CR + 4C) is due to 
chance is less than 0.05. 

On Table II we also give the fraction f for all quasars with known 
redshifts. The difference between f for the optically selected quasars 
and quasars found in radio surveys is similar to the fractions discussed 
above. Since the sample of all known quasars is larger the probability 
that the difference in these samples (dividing the sample at z as 0.4) is 
due to chance is <.0005, a probability less than the probabilities obtained 
for the smaller but complete samples. 

In summary we find that both luminosity fractions of equation (31) and 
(34) encounter difficulties when extended to low radio luminosities. Equa- 
tion (31) disagrees with both the observed redshift distribution and the 
luminosity function of optically selected samples. Equation (34) agrees 
with the observed optical luminosity function but it disagrees with observed 
redshift distribution and predicts higher radio luminosities for quasars in 
optically selected samples than indicated by one set of observations (Wardle 
and Miley 1971). 

We consider these evidences additional support for our earlier conjec- 
ture (Petrosian, 1972, see also the discussion at the end of Section III) 
about the presence of two classes of quasars. 
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VI. Calculation of V/V for Optically Selected Quasars 

m 

In this section we present calculations of ( V/V ) for the Lynds et al 

m 

(1972) sample of quasars which is complete to an infrared magnitude of i = 
16.65. In previous calculations ( V/V ) (Schmidt, 1968, Lynds and Wills, 
1972) the sample was chosen to be complete to a given intensity at a rest 
frame wavelength of 2500 A after corrections for the effects of discrete 
lines. In this procedure few objects are eliminated from a sample that is 
complete to a given apparent magnitude. To avoid this loss of objects, a 
more complicated procedure could be used. We describe this procedure below. 

Let us assume that the spectral luminosity F(v) of an object can 
be described by a power law nonthermal source plus contribution from 
various lines of frequency . Then 

f(v) = F K)[(^r) ° + Z)(^r) ^(v-v.)] , (te) 

f" 

where g(v-V^) describes the line shape ( I g(x)dx=l) and w. is the 

— 00 

equivalent width of the line. If the source is at a redshift z then the 
observed intensity at a frequency v is 

i+z)- v .]J , (43) 

the observed equivalent width is wV(l+z) and the observed magnitude X for 
a filter response S (v) is 

A. 

x = -2.5 log j I(v)s x (v)dv . (44) 

o 

Since the line widths are in general much smaller than the filter 
response widths, we can approximate the line profile by a delta function 
and set v q equal to the center frequency of the filter to obtain 


I(v) = 


F(v q )(1+z) 


l/v( 1+z) \ -a o V''/ V i\ -0 k , 

k \ o / 'o' 


2 6 



icf 0 ^ , F(v, ^ 


— I ('t) ° s x (v)dv [ 1 + Z c i (z >] ' (45) 


WVC 14 - 2 )^ w i s x [ v i /( 1+z ^ 

55 PS — — ————— 

(l+z) 1 “ a oJj^y Q S x (v)dv ( 1 +z) J S x (v)dv 
The last approximation is valid only if filter half power width Av Y « V . 

A O 

For a sample complete to an observed magnitude X the maximum red- 

o 

shift z up to which the objects would be observable is obtained from 
m 

equations identical to equations ( 4-5') and ( 46) with X replaced by X q and 
z replaced by z . We can therefore write 


c t (*) 



\ 

0.4(x o -x) 


21 og 


/z \ 1 +z 

( t ) + ( a o -1 ) log m 


1 + 


1 +z 


+ log 


1 + 


E c i( z ) 

i 

£ C i (z m> 


(Vr) 


Li J 

In general only one strong line contributes to a given filter but different 
lines (if any) may contribute at z and z^. Furthermore, if the observed 
equivalent width ^is much smaller than the filter half power width 

Av y the last term in equation (4-7) can be approximated as 

A 


log(l + C.(z)) 


1 W i S x^ v i / ^ 1+z)) 


2.3 


( 1 +z) J S Y (v)dv 


m 


In order to calculate the values of z^ we need to know the equivalent 

width of the strong lines and the response function S ( v) . Since the sample 

under consideration is limited to infrared magnitude, i.e. X = i, we use 

the response function for this filter published by Braccesi et al. (1970). 

o 

As this response function drops off rapidly for wavelengths \ > 78 OO A 

the only strong lines effecting this filter (up to a redshift of 3 ) are 

o 

Ha, Hp, Mgll and CIII at wavelengths 6562 , 4861, 2798 and I 9 O 9 A respectively. 
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We assume that the equivalent widths are the same for all objects and we 

use equivalent widths (in wavelength units) of 550, 163, 47, and 14 A, 

respectively. These numbers were obtained by straight averaging of all 

rest frame equivalent widths, w^, observed by Oke, Neugebauer, and 

Beklin (1970). Note that for calculating the value of z by this method, 

a trial and error procedure must be used. Once z is obtained we can 

m 

then calculate the values of V/V . 

m 

The results of these calculations are presented on Table III for two 

assumed values of the spectral index a and for the case where we have 

o 

neglected the effect of the lines. As can be seen from this table there 
is a significant difference between the values of ( V/V ) for high and 
low redshift (or luminosity) objects. 

Columns 8 and 9 give the absolute luminosity of the sources calculated 
at two different wavelengths. The numbers in column 7 give the approxi- 
mate infrared luminosities and are obtained from 

l0g infrared = " log ^/(l+z) 1 "^ -0.4(i+Ai)+C ( 49 ) 

where Ai corrects the infrared magnitude for the discrete lines, C is 
an arbitrarily normalization factor and oc =0.5. The numbers in column Q 

O ' 

give the approximate values of logarithmic luminosities obtained from a 
relation similar to equation (49) with the infrared magnitude i replaced 
by the ordinary visual magnitudes, m , with Ai = 0, and with normaliza- 
tion C chosen such that m^ = l8 1 D4 corresponds to log L = 23 w/Hz for 
an object at z = 1 (cf. Schmidt, 1968 ) . These numbers do not exactly 
correspond to log L (where L is the optical luminosity at 2500 A in units 
of w/Hz), their deviation from log L is less than 0.1, an accuracy which 
will suffice for the purpose of our discussion here. The value of the 
luminosities appear to divide the sample into two groups. The low redshift 
( z :£ .4) objects all have log L ~ IQ 22 


except for B272 which has narrow 



emission lines, Lynds et al. 1972, and will be neglected from the remainder 
of our discussion) while the large redshift objects have log L ^ 22.8 ex- 
tending up to log L = 24.0. 

The last two columns on Table III give the inverse of the quantity 



which are needed for calculating the luminosity function (equation ll) . 

c 

The numbers in column 10 are obtained with p(z) = (l+z) for the whole 

sample (cf. Table i) while the numbers in column 11 are calculated assuming 

p(z) = const, for z ^ 0.4 and p(z) = [V(z)/V(0.4)] a for z ^ 0.4. 

Since ( V/V ) for the objects with z > 0.4 is 0.75 (cf . columns 5, 6, 7 

of this table) we have chosen a = 2 (according to equation 18). The ex- 

ponential evolution laws (equation 20) will give numbers for 1/V* very 

m 

similar to the ones in column 11 (cf . our discussion in Section III). 

The cumulative luminosity functions obtained from columns 10 and 11 
are plotted on Figure 6. The important features of the luminosity func- 
tion is the steep rise at high luminosities in contradiction with the flat 

optical luminosity fraction of equation ( 31 ) and a rapid increase at 
22 

L ~ 10 . The latter increase could be evidence for the presence of the 

class I quasars discussed earlier. If so, we could then calculate the 
ratio (pj/pjj.) in equation (20). From Figure (6) we find Pj/pjj ~ 10 . 
Note that at a redshift z ~ 2, Pj/[ ] < 1 which may indicate 

that all of the class I quasars were active class II objects at earlier 
epochs . 
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VII. Discussion 


We have assumed the redshift of quasars to be cosmological and con- 
structed simple laws for the luminosity function and the statistical evolu- 
tion of these objects. We have also assumed density evolution. The 
models developed have come about from the consideration of the difference 
between optically selected quasars (i.e. quasars of all radio luminosities) 
and quasars found in radio surveys which have relatively higher radio 
luminosities. These differences are: 

(1) Dependence of on luminosity. 

(2) Steeper optical luminosity function for optically selected 
quasars. 

(3) Difference between the redshift distributions, in particular 
at low redshifts. 

(4) Smaller radio to optical luminosity ratio for optically 
selected objects. 

Clearly there exist many possible models which could account for these 
differences. Schmidt ( 1970 ) has suggested that the radio and optical lumi- 
nosities are correlated (as in equation 34). This explains the difference 
in item ( 2 ) [and possibly item (4) if the reported (Schmidt 1972) detection 
of radio signals from some of the optically selected quasars is confirmed] 
but cannot account for the other differences. Another possibility is modi- 
fication of the luminosity function of equation ( 31 ) such that the optical 
luminosity function becomes steeper at low values of F (see equation 4o) . 
This accounts for the differences in items (2) and (4). However, the exis- 
tence of a relatively large number of low redshift objects among optically 
selected quasars and the difference between the dependence of ^V/V ) on 
luminosity (if future data show this difference to a more significant level 
than possible at the present time) imply the need for additional modification 
of the luminosity function. 

One possibility suggested here and earlier (Petrosian 1972 ) is that 
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these differences are due to the presence of two classes (I and II) of 

quasars. (Note that this in turn implies that the assumption of pure 

density evolution is incorrect.) The class I quasar typically of 
22 

L ~ 10 w/Hz are more numerous at the present epoch (by a factor of 

23 24 

~100) than the class II quasars of typically L ~ 10 to 10 (Sec- 

26 

tion VI). The class II objects are strong radio emitters F ~ 10 to 
28 

10 w/Hz and show strong density evolution. These are the objects found 

in radio surveys. The class I objects on the other hand are weak radio 
25 8 

emitters f < 10 * w/Hz and show little or no evolution. These are 

the low redshift objects found in optically selected samples. 

According to some definitions of quasars the class I objects may not 
be classified as quasars. They may be bright Seyfert or N-type galaxies. 

They may be objects unrelated to quasistellar radio sources (class II 
quasars). If, on the other hand the spectroscopic and photometric simi- 
larities between the two classes are considered as evidence for some kind 

of relation between the two classes, then the data discussed here suggests 

) 

the following picture. The quasars (class I ) in their quiescent state are 

25 8 22 

low luminosity F = F < 10 ' w/Hz, L = L ~ 10 w/Hz nonthermal 

Q ~ W 

sources. Occasionally there occur explosions in each source with up to 
a hundred fold increase in its optical and even larger increase in its 
radio emission. The result is a class II object. If during these explosions 


L(t) = L ? + 




( 51 ) 


then the statistical correlation between L and F for radio active quasars 
can be explained. Furthermore, if the rate and/or duration of such ex- 
plosions are epoch dependent, then we expect (v/V ) ^ 0.5 for class II 
objects and (v/V ) ~ 0.5 for class I objects. (In fact ( v / v m ) <0.5 if 
these objects could be observed out to a large redshift.) This implies that 
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o 

the ratio of the local densities of the two classes p^/p > 300 to 10 

depending on the density evolution law of class II sources. The observed 

2 

value of this ratio is about 10 . (See Figure 6 and our discussion at the 
end of Section VI.) This value is obtained from extrapolation of densities 
of Class II objects from redshifts of order unity to zero, so that this dis- 
crepancy does not rule out the model. 

If the model is correct than the slope of magnitude counts will be 
< 0.6 and ~ 0.75 f° r class I and class II objects, respectively. Since the 
two classes contribute to the counts approximately equally at 18 magnitude, 
there should be a change in the slope of the optical counts around this 
magnitude. Comparison of source counts below 19 magnitude (Sandage, 1965 ) 
with the counts above this magnitude (Braccesi and Formiggini, 1969 ; 

Sandage and Luyten, 1969 ) indicate the possibility of an increase in the 
slope at higher magnitudes. The presence of such a break in the magnitude 
counts will constitute additional evidence for the model proposed here. Of 
course, data on optically selected samples with larger number of objects 
and on samples limited to lower radio flux densities (S o (l78 MHz) < l) 
will be useful in testing the model . Further radio observations of optically 
selected quasars, and in particular limits on the radio flux density of the 
objects in Schmidt's (1970) sample will also be helpful. 
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Table I. Observed and Calculated Values of <V/V 
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Table II. Fraction of Sources with Redshift Less Than 0.6 
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Corrected for the selection effect due to high v-i values of the low redshift objects. 
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Appendix I 


V/V^ Analysis for Samples with Two Limiting Fluxes 

In the case of two limiting fluxes the relations between the ^V/V^, 

and the general form of luminosity function ¥(L,F,V) (with volume V(z) 

instead of z as epoch indicator) is more complicated. For a given radio 

and optical limit S q and I q , the sources can be divided into two groups, 
F S o 

Those with — > — = R q or the optically limited objects, and the radio 

o „ S 
Y o 

limited objects with — < — . Each group then has a single limit and can 

Li I 

o 

be treated as described in the text . The total number of sources in the 
sample is 


N = N ^ + N , 
opt rad 


(Al) 


where the subscripts refer to the limiting flux. 


N 


opt 


| » opt (L)dL=J dF [ "°dt, r 

A A „ 


„°° J/R_ „v (l/i ) 

!U ° Y(L,F,V)dV , (A2) 


00 


N 


rad 


J oe LR V (F/S ) 

N rad (L)dL = J dLr°dFj m ° Y(L,F,V)dV , (A3) 

o 1 o*o "o 


where V^L/I^) (or V m (F/S Q ) are the maximum values for an optically 
(or radio) limited object with optical (or radio) luminosity L (or F) . 
The average value of V/V is therefore given by 


< dL ir° 7ikf m(F/so) f(L - r ’ v)vdv 


(k) j ° 

f dq J 

LR Q m^I I O 


.V (L/I ) 
nr o' 


Y(L,F,V)VdV > . (A4) 


For pure density evolution 


Y(L,F,V) = Pl (V)Y(L,F) = Pl (V) N ( L ' F ) - , 


V = V (F/I ) , F/L < R 
m m x o' 7 o 


f m 

J 01 

•'o 


(A5) 


(v)dV 


V = V (L/I ) , F/L > R . 
m m v o' ' o 


37 



so that, equation (A. 4) is modified as 


N(L ’ F)x [ v m(l;)] dF + x [ v »(r)] / N < L - F > dF ]- < A6 > 

*o ^*o LR 


where 



as in equation (13) . For discrete sources equation (a 6 ) reduces to 


(*7) 


< V/V m) - 5 t < A8 > 

i=l 

as in equation ( 12 ). 

If x ( v m ) i s independent of V m (dx(V m )/dV m =0) equation (A7) and (13) 

imply 

f/v) - V a , ^ = g§ (A9) 

as in equation (l 8 ). In Section IV we have refered to variation of 
with luminosity. For objects of given optical luminosity 


<kX =[{ LH ° N(L »x[v B (|-)] d r + *[\(£)]£ N( L > )dF ]/|° 


N(L,F)dF. (A10) 


o 

In general, independence of this quantity from L does not imply density 
evolution law of equation (A 9 ) . However, if equation (A 8 ) is valid, then for 
optically limited objects the value of ^7/V^ ^ = x[V m (L/I Q )] becomes 
independent of L. Clearly because of the complete symmetry between L and F 
the same applies to radio limited objects. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 
Figure 5 


Figure 6 


. Comparison of various density evolution laws. For each law the 
exponent is chosen such that p(z=3) = 10. 

. Variation of expected value of V/V with the maximum redshift 

m 

(or luminosity) for various density evolution laws. 

The cumulative luminosity function $(F,L) from the table of 
Lynds and Wills (1972). (a) The radio luminosity function; 
the symbols + ,0,x and stand for log L = 22.6, 23.0, 23.4 and 
22.8, respectively. (b) The optical luminosity function. 


Same as in Figure (3), except for §'(R,L), R = F/L. 


Schematic redshift distributions for samples limited to various 
radio flux density S q . (a) Expected from the luminosity 

function of equation (31)- (b) Expected from the luminosity 


function of equation (34). z = (L /I ) 

, /0 V ' opt v min o' 

z = (R . L . /S ) 1/2 . 
r ' min min o' 


1/2 


z = (F /S ) 
r ' o o' 


1/2 


The cumulative optical luminosity function for optically selected 
samples of Lynds et al . (1972). Dotts from column 11 and crosses 
from column 10 of Table III. 
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Figure 2 
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Figure 





Figure 5 
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